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Abstract 

Continuing previous work, forward-backward multiplicity correlations are studied 
in asymmetric collisions in the framework of the weighted superposition mecha- 
nism of different classes of events. New parameters for the asymmetric clan dis- 
tribution and for the particle leakage from clans in one hemisphere to the opposite 
one are introduced to effectively classify different classes of collisions. This tool 
should be used to explore forward-backward multiplicity correlations in AB and 
pA collisions in present and future experiments at RHIC and LHC. 
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1 Introduction: the weighted superposition mechanism of different classes of events 

It has been found quite recently [[]]] that the weighted superposition of different classes of events 
with negative binomial (NB) properties (events with and without mini-jets in pp collisions, 2-, 
3- or more-jets samples of events in e + e~ annihilation) reproduces, in the GeV energy range, 
the available experimental data on forward-backward (FB) multiplicity correlations of the two 
types of collisions; it is an intriguing result which sheds new light on long range properties in 
multiparticle dynamics and outlines the experimental consistency of the weighted superposition 
mechanism of different classes of events. In addition, the same superposition mechanism has 
been shown to provide in the GeV region an interesting phenomenological tool in order to 
describe the observed shoulder effect in n-charged particle multiplicity distributions (MD's) and 
n-oscillations in the ratios of n-particle factorial moments to the n-particle factorial cumulant 
moments [0, 0, |J. Effects on which QCD has not (up to now) satisfactory predictions 

In the case of pp collisions, essential part of the theoretical background in the GeV region 
has been the assumption that, in each substructure or component described by a negative bino- 
mial (Pascal) multiplicity distribution [NB(P)MD], independently produced clans (they follow 
a Poisson distribution) are binomially distributed in the two hemispheres and that logarithmi- 
cally produced charged particles within each clan distribute themselves again binomially in 
the two hemispheres with an energy independent 'leakage parameter.' This parameter con- 
trols the number of particles, generated by clans in one hemisphere, falling in the opposite one 
and was determined in [[[]] from the data at 63 GeV and 900 GeV. Accordingly, the correct 
reproduction in the GeV region of the experimentally observed increase with cm. energy of 
the forward-backward multiplicity correlation strength as well as of the relation between the 
average charged particle multiplicity of particles lying in the backward hemisphere versus the 
number of charged particles lying in the forward hemisphere (and vice versa) support strongly 
our approach. 

In the case of e + e annihilations, in addition, it should be pointed out that under the same 
assumptions of binomial distributions of clans and of particles generated by clans in the for- 
ward and backward hemispheres, the relatively small value of the forward-backward correlation 
strength in the total charged particle MD and the absence of FB multiplicity correlations in the 
separate 2- and 3-jet samples of events, measured by OPAL collaboration [Q], have been also 
correctly reproduced. Forward-backward correlations in the total charged particle multiplicity 
distribution are indeed only due here to the superposition effect of the two samples of events. 
In our approach, in fact, forward-backward multiplicity correlations in the two separate sample 
of events, each of NB type, turn out to be zero, as no leakage is found from one hemisphere 
to the other in the two separate samples, and in the total sample resulting from the weighted 
superposition of two NB(P)MD we predict a correlation which coincides with the experimental 
one within experimental error [JI[] . 

The striking difference in FB charged particle multiplicity correlation strength for the total 
sample of events between proton proton collisions in the GeV energy domain and e + e~ annihi- 
lation at LEP energy outlines the deep link between FB multiplicity correlations and long range 
correlations, which are expected to be quite strong in the first case and relatively weak in the 
second one. In addition, the observed lack of FB multiplicity correlations in the two separate 
2-jet and 3-jet samples of events in e + e~ suggests that the weak FB correlations seen in this 
reaction in the total sample of events are entirely due to the superposition of the two separate 
samples. This fact can be considered indeed an experimental evidence of the weighted superpo- 
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sition mechanism characteristic of our approach (justified up to now only on the basis of a quite 
successful description of collective effects in high energy energy phenomenology which QCD 
is unable at present to describe) and of the presence or absence of charged particle leakage from 
one hemisphere to the opposite one. 

In terms of clan structure analysis of the NB behaviour for each class of events, it is clear 
that charged particle leakage as well as FB multiplicity and long range correlations are expected 
to be stronger when particle population per clan is larger, a phenomenon which usually goes 
together with the reduction of the average number of clans. 

It should be pointed out that the occurrence of a larger average number of smaller size 
clans and of a smaller average number of larger size clans has a suggestive interpretation at 
parton level, in the framework of a two-step mechanism, and could be related to smaller and 
larger colour flow densities respectively. In a NB description of final charged particle MD for 
a single component (soft or semi-hard in pp collisions, 2- and 3-jet samples of events in e + e~ 
annihilation), after using generalised local parton hadron duality (GLPHD), the first step of the 
parton production process is dominated by the A q ^ q+g (N co i OT , e) vertex and the average number 
of partonic clans corresponds to the average number of bremsstrahlung gluon jets (BGJ), which 
are effective independent intermediate gluon sources (IIGS), and the second step is controlled 
by the vertex A g ^ g+g (N color , e), the gluon self-interaction vertex, whose increase corresponds 
to an enhancement of parton cascading from the IIGS [|J [7p. iV co i or is the number of colours and 
e the fixed cut-off regularization prescription of the theory. 

Coming to the difference between e + e~ and pp reactions, it is likely to expect that the 
observed increase of the average number of clans in e + e~ with respect to the pp case is a con- 
sequence of a stronger activity of the first vertex with respect to the second one, and that the 
opposite will occur for the increase of the average number of partons per partonic clans. It is 
quite clear that the production of a large average number of partons per clan is a consequence of 
a longer cascading process, originated by IIGS generated at relatively high virtuality in regions 
where, being the coupling constant smaller, stronger colour flow between partons should be at 
work (a situation favoured in pp collisions). When the average number of partons per clan is 
relatively small, IIGS are expected to start to be effective at lower virtuality, their cascading 
becomes shorter and colour exchanges reduced with respect to the previous case (a situation 
favoured in e + e~ annihilation). It seems therefore that the occurrence of stronger FB multi- 
plicity correlations, long range correlations and related particle leakage enhancement from one 
hemisphere to the opposite one at hadron level are a specular image of stronger cascading from 
high virtuality IIGS and of larger colour exchange in this region at parton level. 

We believe in fact that the understanding of FB multiplicity correlations at hadron level is 
a possible starting point in order to study new effects of colour quantum number exchanges in 
multiparticle dynamics. 

Accordingly, we decided to continue our search initiated in Ref. [jlj], where FB correlations 
have been understood in the GeV region in symmetric reactions (like pp collisions and e + e~ 
annihilation) and for symmetric definition of the hemispheres by assuming at hadron level: 

a. NB behaviour for the (forward plus backward) MD of each component or substructure 
(class of events), i.e., clan structure analysis is assumed to be applicable to each of them. 
The generalisation to the class of compound Poisson distributions (CPD) is of course 
possible: the generating function Gqpd(z) in this case can be written as follows 

G?cFD(«) = exp{JV fl [(/ c (z)-l]} (1) 
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where N g is the average number of generalised clans and g c (z) the generating function 
of the MD of charged particles originated by a generalised clan (e.g., it is a logarithmic 
distribution when the MD for a single component is a NB(P)MD). 

b. clans are therefore independently emitted and their distribution is Poissonian. In case of 
a generic CPD we talk, as previously stated, of generalised clans. 

c. clans are binomially distributed in the forward and backward hemispheres (as are gener- 
alised clans). 

d. logarithmically produced charged particles in each clan are also binomially distributed in 
the forward and backward hemispheres but with different probabilities p and q(p + q = 1, 
p different from q); p controls the leakage from one hemisphere to the other: p — 1 means 
that no particle leaks, while p larger or equal to 0.5 and smaller than 1 indicates leakage. 
In case of a generic CPD, the n-charged particle MD is generated by g c (z), but p and q 
retain their meaning as described. 

In order to perform our calculations in pp collisions one extra assumption has been added, 
i.e., that the particle leakage parameter is constant throughout the GeV region for the two sepa- 
rate classes of events. 

Possible scenarios in the TeV region based on extrapolations from data in the GeV region 
and on the weighted superposition mechanism of soft (without mini-jets) and semi-hard (with 
mini-jets) events have been indeed studied and predictions given on n-charged particle multi- 
plicity distributions and on the ratios of n-particle factorial moments to the n-particle factorial 
cumulants moments general properties in the TeV region . 

An attempt to predict the energy dependence of the forward-backward multiplicity correla- 
tion strength as well as ub^f) vs rip general trends in the new energy range available at CERN 
with Alice detector has been considered [[[J]. The word is now to experiment which is supposed 
to test all these three sets of predictions at 14 TeV. 

The problem we want to face in this paper, on the theory side, is the generalisation of 
the approach discussed in Ref. [|T]] for determining forward-backward multiplicity correlations 
properties for a single component, in more complex asymmetric reactions (like heavy ion AB 
and pA collisions) and to provide a general framework for the study of forward-backward mul- 
tiplicity correlations which includes symmetric reactions (like pp and AA collisions and e + e~ 
annihilation) as a particular case. 

Corner stones of our argument remain of course assumptions a, b, c, d. 

They could hardly be abandoned in view of their success in giving a good phenomeno- 
logical description of available experimental data on FB multiplicity correlations in symmetric 
reactions. 

The generalisation of the approach to asymmetric reactions will concern therefore mainly 
how to implement, in the framework defined by assumptions a.b.c.d, the asymmetry of the 
reaction and the asymmetric definition of the forward and backward hemispheres. 

In view of the lack of sound experimental data on asymmetric reactions and of the related 
analyses in terms of two or more component substructures both in full phase space and in 
rapidity intervals, the present paper should be considered as a stimulus to experimentalists of 
the new generation machines for a deeper analysis of total n-charged particles MD and a more 
satisfactory understanding of forward-backward multiplicity correlations as the cm. energy 
increases and at fixed cm. energy within rapidity intervals. 
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2 The general asymmetric case 

Generalisation of the study performed in Ref. [[]]] on FB multiplicity correlations consists at the 
present stage of investigation in assuming that 

A. the particle leakage percentage from the backward (B) to the forward (F) hemisphere, 
q B , is different from the particle leakage percentage from F to B hemisphere, q F , with 
Pf + qF = 1 and p B + qs = 1, Pf and p B being the corresponding percentages of 
particles not leaking from one hemisphere to the opposite one and remaining in the F and 
B hemisphere respectively; notice that in Ref. [|T]j, Pf = Pb = P- 

B. Poissonianly generated clans are asymmetrically (but binomially) distributed in the two 
hemispheres with the asymmetry parameter r different from 1/2, s being its complement 
to 1 in the opposite hemisphere; in Ref. [[]]], notice that r = s = 1/2 in addition to 

Pf =Pb = P- 

2.1 The generating function 

In general, the joint MD P(n F , »b) for n F particles in the forward hemisphere and n B parti- 
cles in the backward one is related to the global MD P{n) though the probability distribution 
f(n F \n), which gives the probability to have n F F-particles when the total number of particles 
is n, as follows: 

P(n F , n B ) = P(n F + n B )f(n F \n F + n B ). (2) 
The generating function G(z F , z B ) for the joint distribution then satisfies 

G(z F ,z B ) = z n /z n B B P{n F ,n B ) = J2z%P(n)g f (z F /z B ;n), (3) 

riF,nB n 

where gj(z; n) is the generating function for f(n F \n); in case it is the binomial distribution 
with parameter p, 

f(n F \n) = ^p n ni~p) n - nF , (4) 
then its GF is (defining q = 1 — p): 

g f (z;n) = (q + pz) n , (5) 
and one obtains a considerable simplification of Eq. (0): 

G(z F , z B ) = g{z F p + z B q), (6) 
where g(z) is the generating function of the global distribution P{n): 

g(z)=J2z n P(n). (7) 

n 

Accordingly, we proceed now to calculate the GF for the general case. The formulae below 
are heavily based on previous work [p. We consider a clan with logarithmic MD produced in 
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the F hemisphere and assume that each particle has the same probability p F not to leak into 
the B hemisphere. Then the distribution at fixed number of particles is binomial, and we have 
immediately an application of the just-explained scheme: the GF for the joint distribution within 
a F-clan is thus 

g c ,F(z F , z B ) = g\ og {z F p F + zbQf) i (8) 
for a clan produced in the B hemisphere we have the corresponding GF: 

9c,b(z f , z b ) = 9\ og {z F qB + zbPb)- (9) 
Here the GF for the logarithmic distribution of parameter (3 is 

^ g (*)=log(l-z/3)/log(l-/3), (10) 
where, in terms of standard NB parameters n and k, 

n + k 

and is related to the average number of particles per clan via the formula: 

nc = (p- 1)4(1-/3)- (i2) 

Recalling now Eq. (34) of [p, it is easy to calculate the GF of the joint distribution of F and 
B particles at given numbers N F of F clans and N B of B clans: 

g(z F , z B \N F , N B ) = [g c , F {z F , z b )} Nf [g c ,B{z F , z b )] Nb : (13) 

due to the fact that all clans are by definition independent from each other, we can convolute 
the respective MD's, which corresponds to multiply together the GF's. 

In order to sum over the clan MD, let us remember that we have assumed that clans are 
Poisson distributed and independent of each other, thus the joint distribution can be written as: 



where N = N F + Nb is the average number of clans, and r = N F /N is the fraction of clans 
emitted in the F hemisphere. The corresponding GF is 

g{z F ,z B ) =exp{N[rz F + (l-r)z B -l]} . (15) 

We can now perform the last step in the calculation, again exploiting the general properties 
of the binomial distribution: 

g(z F ,z B ) = ^2^2g(z F , z B \N F , Nb)V(N f , N B ) = Q (g c , F (z F , z B ), g c ,B(z F , z B )) ■ (16) 

Np N B 

g(z F , z B ) = exp {rN[g Xog (z F p F + z B q F ) - 1]} exp {(1 - r)N[g Xog (z F q B + z b Pb) - 1]} ■ 

(17) 

Notice that the above formula is valid for the NB(P)MD, but can easily be extended, as antic- 
ipated, to any compound Poisson distribution, provided all correlations are exhausted within a 
(generalised) clan: it is sufficient to replace g\ og with the appropriate GF within a clan, g c (z). 
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2.2 The correlation strength and the marginal distributions 

We can now use the above result to deduce the correlation strength b: 

_ (n F n B ) - n F n B 

~ [((n F )-n F )((nl)-nl)} 1/2, 

since everything can be read from the GF: 

d 2 g 



(n F n B ) 



dz F dz B 



(18) 



(19) 



Zp=Z B = l 



71: 



dg 

dzi 



(n 2 ) 



d 2 g 



Zp=Z B =\ 



dzf 



(20) 



Z F =Z B = 1 



with i = F, B. 

The average number of F-particles at fixed number of B-particles is also easily obtained 
though differentiation, since the MD of n F at fixed n B is given by 



p{n F \n B ) 



p(n F ,n B ) 



p{n F ,n B ) 



(n B ) E„ f pK^b)' 
where pinJg(n B ) is the marginal distribution in the B hemisphere; the GF is 

g(z F \n B ) = y^ j z F F p(n F \n B ) 



(21) 



q^9{Zf, zb> 



Qn 



np 



dz n B * 



(22) 



z B =0 



( Ft} 

being g(l, z B ) the GF of phwg. The first moment is obtained from the GF by differentiating 

Qn B q 



once: 



n F {n B ) 



dg(z F \n B ) 



dzi 



dzl B dz F 



g(z F ,z B ) 



Zp = l 



Qn B 

-Qj&g{z F ,z B ) 



(23) 



z B =0, zp=l 



Analogous formulae for the B hemisphere quantities can be easily obtained in the same way. 
The corresponding marginal distributions are g(z B = z, 1) and g(l, z F = z): 

g(z, 1) = exp [rN[g Xog (zp F + q F ) - 1]} exp {sN[g Xog (zq B + p B ) - 1]} ; (24) 

g(l, z) can be obtained from Eq. ( |24l ) by interchanging parameter pi with qi (i=F,B)- 

The marginal distribution of Eq. ( |2"3| ) is the product of the GF's of two NB(P)MD's with 
characteristic NB parameters {nrp F , kr} and {nsq B , ks} respectively as can be seen immedi- 
ately by noticing that 



g(z,i) 



kr 



-kr 



(25) 
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g(l, z) can be defined in an analogous way. 

g(z, 1) and g(l, z) are the GF's of the MD's obtained by convoluting the MD's for particles 
generated by F (B) clans which stay in the F (B) hemisphere and do not leak in the opposite 
hemisphere with the MD's for particles generated by B (F) clans which are leaking in the F (B) 
hemisphere. 

Although g(z, 1) (and g(l, z)) is the product of the GF's of NB(P)MD's it is not the GF 
of a NB(P)MD. This consideration notwithstanding, it is interesting to remark that g(z, 1) (and 
g(l, z)) is an infinitely divisible distribution (IDD), a fact which allows to define the generalised 
clan concept, as already remarked. 

In fact Eq. ( p4| ) can be rewritten as follows 

g(z, 1) = exp {N[rg log (zp F + q F ) + sg Xog {zq B + p B )} - N} (26) 

i.e., as 

g(z, 1) = exp [N(A(z) - 1)] (27) 

with A(z) = [rg\ og (zp F + q F ) + sg\ og (zqB +Pb)) and A(0) different from zero. This last remark 
implies that the probability of generating zero particles is different from e~ N . Therefore the 
probability of generating empty clans is also different from zero, a result in contrast with the 
standard clan definition (each clan contains at least one particle), definition which we would 
like of course to enforce. 

In order to do that, let us add and subtract in the exponent of Eq. ( ]2~7| ) the term A(0), i.e., we 
rewrite Eq. ( j27"l ) as follows 

g(z, 1) = exp {N[A(z) - A(0) + A(0) - 1]} . (28) 

From Eq. ( p8| ) one obtains the compound Poisson distribution belonging to the class of IDD 

g(z,l)=e X p{N g [G g (z)-l]}, (29) 

with N g = [l- A(0)]N and G g (z) = [A(z) - A(0)][1 - A(0)}-\ 

Accordingly, one can define the average number of generalised clans, N g , in terms of the 
standard NB parameters n and k, and of p F and ps'- 

N g = N[l - A(0)] = N{1 - [rg log (q F ) + sg log (p B )]} 

= — In (7(0, 1) = rk\n(n + k — nq F ) + sk\\i{n + k — nps) — kink. 

The backward marginal multiplicity distribution and related properties can be be obtained from 
the forward one by interchanging parameter pi with q { (i=F,B). 

3 The cases of partially removed symmetry 

3.1 Symmetry for clans only (r = s = 1/2, p F 7^ ps, qF 7^ Qb) 

The symmetry of the reaction can be broken partially by assuming that binomially distributed 
clans go fifty per cent in the F hemisphere and fifty per cent in the B hemisphere (r = s = 1/2) 
but particle leakage from clans in F to B hemisphere, q F , and particle leakage from clans in B 
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to F hemisphere, q B , (as their complements p B and p F with p F + q F = 1 and p B + g B = 1) are 
different. In other words, assumption B from the previous section is valid but A is not. 

The forward-backward joint charged particle multiplicity GF for one single component turns 
out to be 



g(z B , z F ) = exp {\N[g Xog (z F p F + z B q F ) - 1]} exp {\N[g log (z B p B + z F q B 
and for the corresponding marginal B distribution (z B = z, z F = 1) 

g(z, 1) = exp {lN[g\ og (PF + zq F ) - 1]} exp {\N[g {og (zp B + q B ) - 1]} . 
In terms of standard NB parameters Eq. ( |32| ) becomes 



1]} 



(3D 
(32) 



1 + 



1 + 



— — V 

1 + t(! - «) + -j^PF(lB> 



1-z) 



-fc/2 



-fc/2 



(33) 



The symmetric definition of the two hemispheres is removed here in one component by as- 
suming that particle leakage parameters in the forward and backward hemispheres are different. 

The GF of the forward marginal multiplicity distribution thanks to the quadratic term as al- 
ready shown in the general case is no more a NB(P)MD, although it is the product of the GF's of 
two NB(P)MD's with characteristic parameters {\np F , and {\nq B , with asymmetric 
average charged particle multiplicities. The GF of the marginal distribution becomes the GF 
of a NB(P)MD in some special cases, i.e., for q B = 0, corresponding to no particle leakage 
from the backward to the forward hemisphere, and for n < fe, a situation which occurs in the 
Poissonian limit for average charged particle multiplicity much less than the k parameter and 
almost coinciding with the average number of clans. The GF g(z, 1), although not NB, is still 
an infinitely divisible distribution, as expected. 



3.2 Symmetry for particles within clans only (r ^ s ^ 1/2, p F = p B , q F = q B ) 

Another way to remove partially the symmetry is to use assumption A without B, i.e., to as- 
sume that binomially distributed clans are not symmetrically subdivided between F and B hemi- 
spheres, but particle leakage from clans in one hemisphere to the other is the same (r ^ s ^ 1/2 
with r + s = 1, and p F = p B = p, q F = q B = q with p + q = 1, and p larger or equal than 1/2 
and smaller than 1). 

The FB joint particle multiplicity distribution GF for one component, g(z F , z B ), becomes 
in this case 

g(z F , z B ) = exp [Nr[g Xog (z F p + z B q) - 1]} exp {Ns[gi og (z B p + z F q) - 1]} , (34) 
and the corresponding forward marginal charged particle MD GF 

g(z, 1) = exp [Nr\g Xog (zp + q) - 1]} exp {Ns[g log (zq + p) - 1]} , (35) 

which in terms of standard NB parameters of the total charged particle MD of the component 
under investigation is 

-kr r inn -I — fc'S 



1 + > 



(36) 
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The forward marginal MD GF, although again not of NB type, is the product of the GF's of 
two NB(P)MD's with parameters {nrp, kr} and {nqs, ks} respectively. The two sets of param- 
eters characterise the group of non-leaking particles lying in the F hemisphere and those which 
leak from the backward hemisphere to the forward one. The backward marginal distribution 
can be easily obtained from the above equations by interchanging parameter r with s. 



4 The symmetric case (r = s = 1/2, p F — p B , q F = q B ) 



In this case, assumptions A and B are both rejected. The formula studied in [JT|] for the joint 
charged particle MD GF in the symmetric case follows from Eq. ( |T7| ) by taking p F = Pb = P, 
Qf = Qb = Q with p different from q, p + q — 1, and r — s — 1/2. 
We get 



g(z B , z F ) = exp {\N[g log (z F p + z B q) - l}[gi og (z B p + z F q) - 1]} 
[k + n(l- pz B - qz F )] k/2 
[k + n{l — pz F — qz B )] k / 2 

and for the corresponding forward marginal distribution g(z, 1) 



[l + f(l--)]- fc/2 [l + f(l 



-ft/2 



n 



1 + -(1 - z) +pq 

fv 



-fc/2 



(37) 



(38) 



In conclusion, Eq. ( |P7[ ) can be considered the wanted generalisation of Eq. (|37|) to the case of 
asymmetric reactions and an asymmetric definition of the forward and backward hemispheres. 



5 Behaviours of n F (n B ) vs. n B and of n B {n F ) vs. n F 

In this section we examine the relation between the average number of particles in one hemi- 
sphere, npins) or fiB(n F ), at fixed value of the number of particles, respectively % or n F , in 
the other hemisphere, according to Eq. (Q). 

Let us start by recalling that the von Bahr-Ekspong theorem [§] implies, since we assume 
the total MD to be of NB type, that there is no linearity in fiB(n F ) vs n F (and in npins) vs 
ub) unless the MD for F-particles (and for B-particles!) at fixed total number of particles is 
binomial, in which case the GF is simply: 

exp[g log (z F p + z B q) - 1]. (39) 

Notice that it does not make sense to distinguish p F from p B in this case; however, if one 
did not distinguish F-clans from B-clans, i.e., if each clan emitted the same fraction p in one 
hemisphere, then in Eqs. (|8|) and @ above one would put p F = q B = p and Eq. ( |T7| ) would 
reduce to Eq. (|39|). 

In the following the symmetry in the clan distribution is contrasted with the asymmetric 
clan distribution. In addition, leakage parameters for the two components are taken either equal 
or different. The resulting pictures of forward-backward multiplicity correlations are shown 
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Figure 1: Behaviours of npins) vs. ub (dashed lines) and of nsinp) vs. n F (solid lines) for r = s = 1/2 and 
values of the leakage parameters as indicated. In order to illustrate a practical case, we have taken for the NB 
parameters n — 39.5 and k — 7.0, corresponding to the soft component of pp collisions at 14 TeV in the scenarios 
examined in Ref. [[||. In panes where only the solid line is visible, it means the dashed line coincides with it. 

in Fig.s [1] and Q of course our choice of the involved parameters is arbitrary. The word is 
again to experiment on symmetric collisions (r = s = 1/2) in order to understand eventual 
deviations from identity of particle leakage parameters from clans in the forward and backward 
hemispheres and to experiment on asymmetric collisions (r ^ s ^ 1/2) in order to measure the 
different leakage parameters in the two hemispheres. All together, expected different values of 
r and s as well as p F and p B parameters could lead to a new intriguing classification, in terms 
of FB multiplicity correlations, of high energy collisions and their substructures. 

5. 1 The case r = s = 1/2 ( symmetry in the clan distribution ) 

We start by examining the case in which r = 1/2, illustrated in Figure [I]. When p F = p B = 1 /2, 
there is perfect linearity in n B {n F ) versus n F . This is in agreement with the mentioned theorem, 
because here we are saying that, within each clan, particles are binomially distributed in F and 
B with the same probability 1/2, thus the fact that particles are produced in clans does not make 
any difference. When 1/2 < p F = Pb < 1, on the contrary, the fact that the production happens 
in two steps becomes again important and the relation between n B {n F ) and n F is non-linear. 
Linearity is recovered again in the limiting case of no correlations (p F = p B = 1), see below. 
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When p F 7^ p B , the symmetry is lost because the particles leaking from the F hemisphere 
into B are not compensated (if p F < p B ), or are over-compensated (if p F > p B ), by those 
leaking from B into F. Thus one obviously finds that n B (n F ) vs n F is not the same as n F (n B ) 
vs n B . Indeed the two marginal distributions differ from each other, but can be obtained one 
from the other by exchanging p F with p B . In general for p F > p B there is less leakage from F 
than from B: particles prefer to stay in the F hemisphere thus one has the line of n F (n B ) vs n B 
above n B (n F ) vs n F . Vice versa, the opposite is true for p F < p B . 

One can notice that when p F = 1, n F [n B ) vs n B is a straight line, and when p B = 1, 
n B {n F ) vs n F is linear. The fact is that when p F = 1, clans that fall in the F hemisphere 
do not contribute to the B hemisphere. The only relation between n B and n F is given by B 
clans leaking to the F hemisphere. And this relation is by definition binomial for each clan: the 
number of clans is Poisson distributed, thus the overall is still binomial with the same parameter, 
as was shown above. Indeed, the GF with p F = 1 is: 

g(z F , z B ) = exp {rN[g log (z F ) - 1]} exp {(1 - r)N[g log (z F q B + z B p B ) - 1]} ; (40) 

applying Eq. (|23|), the first term gives just the constant rn, while the second one is of the same 
type of Eq. ( P^| ) and thus also gives a linear behaviour. Notice that this happens independently of 
the value of r. A corresponding result can be obtained for p B = 1 by appropriately exchanging 
the roles of F and B. 

5.2 The case r ^ s ^ 1/2 (asymmetry in the clan distribution) 

When r increases over 1/2 (e.g., in Figure ^ the case r = 3/4 is illustrated), the excess of clans 
in the F hemisphere implies an increase of n F {n B ) vs n B and a decrease of n B (n F ) vs n F , 
enlarging the differences between the curves when p F > p B and reducing them when p F < p B 
(n F (n B ) vs n B can even become larger than n B (n F ) vs n F , depending on the values of r, p F 
and p B ). Exchanging p F with p B does not restore the symmetry. If p F = p B = 1/2, then the 
value of r has no influence on the outcome: if particles within each clan have a 50% chance 
of going into the other hemisphere, it makes no difference if clans prefer one hemisphere over 
the other, and the two curves coincide again. This is not true if p F = p B > 1/2. Finally, in 
agreement with the remark at the end of the last subsection, we notice that even in the case 
r > 1/2 when p F = 1, n F (n B ) vs n B is a straight line, and when p B = 1, n B (n F ) vs n F is 
linear. 



5.3 Asymmetry in the average number of particles per clan 

All the above has been studied in the framework of a NB(P)MD describing in one component 
the forward plus backward MD. However, there are low energy data on pA collisions in which 
NB behaviour has been found to hold separately in the forward and in the backward hemisphere, 
with different NB parameters [|H3|]. Because clans are independently emitted, the only correction 



to make to our formulae is to allow for different average numbers of particles per clan in the 
two hemispheres: in other words, the parameter [3 of the logarithmic distributions in Eq.s (|]) 
and @ will be different, say (3 F and (3 B . Now, we immediately get 

g(z F , z B ) = exp {rN[g log (z F p F + z B q F ; (3 F ) - 1] } exp {sN[g log (z F q B + z B p B ; (3 B ) - 1] } 

(41) 



12 



120 

90 
60 
30 


120 



p F =0.50 p B =0.50 



120 




_ I I I I I I I I 


| 1 1 1 _ 


: p F =0.50 p B =1.00 








* T ' TtTtTtTttI 


111 = 



20 



40 



60 



80 



p F =0.75 p B =1.00 




20 40 



60 



-= 
80 



60 



3 
80 



_ I I I I I I I I 


l l l _ 


; p F =1.00 p B =1.00 








- i i i I i i i I i i 


ill- 



60 



80 



Figure 2: Same as Figure [l] but for r = 3/4, s = 1/4. 



in place of Eq. (|T7|). The product of these two NB(P)GF is not a NB(P)MD, but it is still, by 
construction, an infinitely divisible distribution. 



6 Conclusions 



It has been shown that assuming different particle leakage percentages (ps ^ Pf) for binomially 
generated particles from clans in one hemisphere to the opposite one and asymmetric (r ^ s) 
distribution in the two hemispheres of binomially generated clans, a general formula for the 
generating function of the joint (rip, n B )-charged particle multiplicity distribution for each class 
of events (or substructure) can be obtained when the total MD GF is of NB type. The formula 
reduces to that discussed in Ref. [jH] for p B = p F and r = s = 1/2. Of particular interest are 
also the cases in which the symmetry is only partially removed (assumption A without B and 
B without A). All above-mentioned results, although explicitly derived for substructures of NB 
type, can be easily generalised to any discrete infinitely divisible MD. This search is relevant 
for the study of forward-backward multiplicity correlations in non-identical heavy ion and in 
proton-nucleus collisions. Accordingly, the newly introduced particle leakage and asymmetry 
parameters can be considered as effective indeces classifying different classes of collisions. 
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